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Fundamental  Studies  in  Crack  Initiation 


ABSTRACT 

Mechanistic  investigations  of  damage  evoiution  before  crack  initiation  in  an 
amorphous  polymer  show  that  damage  consists  of  a  core  of  highly  dense  crazing  and 
a  peripheral  less  dense  zone  of  crazing.  Damage  characterization  is  carried  out  at 
consecutive  configurations  of  the  damage  zone.  Analysis  of  the  kinematics  of  damage 
at  different  times  involves  comparisons  of  the  inertia  moments  of  damage  distributions. 
The  results  indicate  that  damage  evolution  between  consecutive  configurations  can 
be  approximated  by  a  linear  transformation  of  the  space  variables.  Thus,  the  process 
of  damage  growth  can  be  described  by  translation  and  deformation  of  the  damage 
zone.  The  growth  rates  of  the  damage  zone  movements  decrease  until  crack 
initiation.  In  all  cases,  the  average  densities  exhibit  a  damping  type  behavior  with  the 
number  of  cycles.  The  crack  initiates  within  a  core  zone  immediately  ahead  of  the 
stress  concentrator.  The  experimental  results  suggest  that  damage  density  within  the 
core  zone  is  independent  of  the  loading  conditions  considered  herein.  This  value  is 
approximately  equal  to  the  damage  density  around  the  crack  tip  during  slow  crack 
propagation.  The  crack  length  at  initiation  is  found  to  increase  exponentially  with  the 
stress  level.  A  simple  decaying  exponential  relationship  relates  the  crack  initiation 
times  and  the  applied  stress  level.  This  result  Is  consistent  with  the  fracture  models 
based  on  absolute  reaction  theories. 

An  energy  release  rate  due  to  damage  growth  before  crack  initiation  is  evaluated  with 
the  use  of  a  semi-empirical  method  for  energy  calculations  and  experimental 
measurements.  Correlation  of  the  elementary  movements  of  the  damage  zone  with 
the  energy  release  rate  shows  that  the  rate  of  damage  growth  decreases 
monotonically.  This  is  consistent  with  the  growth  behavior  of  an  average  damage 
density  within  the  zone.  The  experimental  data  have  shown  that  crack  initiation 
occurred  when  damage  density  within  a  core  zone  ahead  of  the  notch  tip  reached  a 
critical  level,  independent  of  loading  conditions.  Assuming  that  damage  evolution  is  a 
stress  -  temperature  driven  process,  it  is  shown  that  a  first  order  reaction  equation 
describes  reasonably  well  damage  growth  within  the  core  zone.  The  rate  constants 
are  evaluated  by  invoking  principles  of  kinetic  theories  of  fracture  and  experimental 
results  on  crack  initiation  times.  An  activation  energy  for  defect  nucleation  is  found 
half  the  energy  for  thermal  destruction  and  twice  the  enthalpy  of  activation  for 
secondary  chain  motions  below  glass  transition  temperature.  This  result  Indicates  that 
chain  motions  and  chain  scission  are  possible  rate  determining  process  leading  to 
crack  initiation.  This  is  consistent  with  the  observation  that  dense  crazing  precedes 
crack  initiation. 
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SIGNIFICANT  FINDINGS 


1 .  A  damage  zone  accompanied  the  process  of  crack  initiation.  Within  the  zone 
two  patterns  of  damage  were  distinguished:  a  core  of  highly  dense  damage  and  a 
peripheral  less  dense  crazes  around  the  core  zone.  Crack  initiation  occurred  within 
the  core  zone. 

2.  Damage  evolution  before  crack  initiation  evolved  in  a  self  -  similar  manner.  That 
is,  a  linear  transformation  related  the  points  of  equal  damage  at  consecutive 
configurations  of  damage. 

3.  While  the  growth  behavior  of  the  damage  zone  center  depended  upon  the 
loading  history,  the  average  rates  of  the  damage  zone  expansion  and  distortion 
decreased  until  crack  initiation.  Moreover,  the  dependency  of  an  average  damage 
density  with  the  cycle  number  exhibited  a  damping  behavior. 

4.  Damage  density  within  the  core  of  damage  was  found  independent  of  the 
loading  conditions,  however,  the  pattern  of  the  peripheral  crazes  was  dependent  upon 
the  loading  conditions.  The  large  difference  in  crack  initiation  times  was  attributed  to 
the  role  of  time  dependent  processes  within  the  damage  zone  and  the  interaction 
between  the  crazes. 

5.  The  crack  length  at  initiation  was  found  to  increase  exponentially  with  the 
applied  stress. 

6.  The  time  to  crack  initiation  was  related  to  the  stress  level  with  an  exponentially 
decaying  relationship.  The  results  agree  with  fracture  models  based  on  reaction 
theories. 

7.  Using  a  number  of  assumptions  regarding  craze  openings  and  a  semi  •  empirical 
method,  an  energy  release  rate  corresponding  to  damage  growth  before  crack 
initiation  was  evaluated.  Correlation  of  the  energy  release  rate  with  the  elementary 
movements  of  the  damage  zone  showed  that  the  overall  behavior  of  damage  growth 
before  crack  initiation  decreased  with  energy  release  rate.  This  trend  is  consistent 
with  the  growth  behavior  of  an  average  damage  density  before  crack  initiation. 

8 .  The  data  showed  that  damage  evolution  within  the  core  zone  can  be  approximated 
by  a  first  order  kinetic  equation. 

9 .  An  activation  energy  for  defect  nucleation  was  found  half  the  activation  energy  for 
thermal  destruction  and  about  twice  the  enthalpy  of  activation  for  secondary  polymer 
chain  movements.  This  result  suggested  that  secondary  chain  motion  and  scission 
are  the  rate  determining  processes  leading  to  crack  initiation. 

The  results  of  the  present  studies  and  those  reported  in  [17,  18]  indicate  certain 
similarities  between  the  characteristics  of  damage  before  crack  initiation  and  during 
slow  crack  growth.  These  are,  (a)  damage  evolution  can  be  described  by  a  linear 
transformation,  (b)  the  densities  of  damage  within  a  core  zone  at  crack  initiation  and 
during  slow  crack  growth  are  of  the  same  level. 
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1.  INTRODUCTION 

Prediction  of  life  time  of  a  structural  component  subjected  to  fatigue  or  creep  loads 
assumes  a  position  of  prominence  in  engineering  design.  It  becomes  even  more 
important  when  considering  the  continuous  increase  in  complexity  of  engineering 
structures  and  the  large  costs  resulting  from  premature  failures.  Thus,  characterization 
of  the  mechanisms  leading  to  crack  initiation  and  growth  is  an  issue  of  central 
importance  to  engineers  and  material  scientists. 

Contemporary  approaches  to  life  time  prediction  consider  a  fracture  process  in 
three  phases:  crack  initiation,  slow  crack  growth,  and  fast  fracture.  The  latter  phase  is 
very  short  in  duration  and  thus,  the  time  of  crack  initiation  and  that  of  slow  crack  growth 
account  for  the  useful  time  of  a  structure. 

Significant  efforts  have  been  devoted  to  studying  slow  crack  propagation  in 
various  materials.  This  is  demonstrated  by  the  vast  amount  of  experimental  data  and 

several  models  that  have  been  proposed  to  describe  crack  growth  with  the  stress 
intensity  factor  K.,,  or  energy  release  rate  [1,  2].  These  analyses  play  an  important 

role  in  life  time  prediction.  However,  experimental  and  theoretical  studies  of  crack 
initiation  have  received  limited  attention  in  the  literature.  Moreover,  fracture 
mechanics  approaches  to  fatigue  do  not  consider  crack  initiation  [1].  Instead,  the 
concept  of  crack  growth  threshold  is  employed  [3,  4].  Values  of  fatigue  threshold 
Indicate  the  cyclic  stress  intensity  factor  AK^f^,  below  which  long  cracks  remain 

dormant.  The  corresponding  AK^|^  is  taken  as  a  material  parameter.  Although  the 

study  of  fatigue  thresholds  and  near  threshold  behavior  has  highlighted  the  effects  of 
load  history  in  influencing  crack  growth,  and  has  provided  insight  into  the  different 
mechanisms  of  crack  closure,  the  threshold  idea  has  met  with  limited  success  in 
engineering  design.  This  stems  from  the  fact  that  the  use  of  a  fatigue  threshold 
represents  a  conservative  design  criterion  and  from  questions  regarding  as  a 

material  parameter  [5-7]. 

Mechanistic  investigations  in  various  materials  have  shown  that  fatigue  crack 
initiation  and  growth  are  preceded  by  various  forms  of  damage  that  nucleate  at 
heterogeneities  within  the  material.  In  metals  initiation  is  related  to  intense  slip 
processes,  extrusions-intrusions  and  persistent  slip  bands  [8,  9].  Microcracks 
nucleated  at  inclusions  or  second  phase  particles  in  the  case  of  smooth  specimens 
have  also  been  reported  by  several  researchers.  Observations  of  the  same  process  in 
both  notched  and  smooth  specimens  of  polymers  have  emphasized  the  role  of  crazing 
in  amorphous  and  semi  -  crystalline  polymers  [10-13].  Damage  nucleation  and 
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growth,  whether  in  the  form  of  microcracks,  voids,  homogeneous  transformations, 
crazes,  etc.,  are  processes  that  absorb  energy,  energy  that  otherwise  would  be 
available  to  drive  the  crack.  Therefore  the  nature  and  extent  of  damage  determine  the 
time  to  crack  growth  and  toughness  of  the  material.  Under  these  circumstances,  use  of 

fracture  mechanics  parameters  may  not  be  always  justified.  That  is,  when  the 
restrictions  of  'K-dominance'  or  'J., -dominance'  are  not  met,  K.,  and  fail  to  correlate 

with  the  fracture  process.  Examples  are  large  scale  plasticity,  microcrack  interaction, 
the  behavior  of  small  cracks,  etc. 

Recognizing  the  importance  of  damage  on  crack  initiation  and  growth,  several 
analytical  approaches  have  been  developed  to  account  for  its  effects  on  stress 
intensity  factor  and  energy  release  rates  (for  a  recent  review  see  [1 4]  and  references 
therein).  It  is  difficult,  however,  to  assess  the  contributions  of  crack  damage 
interactions,  describe  the  local  kinetics,  and  the  effects  of  the  microstructure  when  a 
large  number  of  microdefects  is  present.  These  powerful  mechanisms  are  the  main 
obstacles  in  modeling  fracture  in  a  large  class  of  modern  engineering  materials. 

One  way  to  treat  fracture  phenomena  where  dense  damage  accompanies  crack 
initiation  and  growth  is  to  use  the  ideas  of  self  -  similarity.  Drawing  on  the  success  of 
statistical  self  -  similarity  of  vortex  cascade  in  understanding  turbulence,  it  is  proposed 
that  the  same  principles  can  be  applied  to  damage  and  fracture  of  engineering 
materials  [15,16].  Indeed,  experimental  observations  in  a  number  of  materials  have 
shown  that  damage  growth  displays  certain  invariant  characteristics.  For  example,  the 
evolution  of  a  normalized  pore  size  distribution  in  polycrystalline  steel  under  creep 
follows  a  pattern  that  is  independent  of  time  and  strain  level  [15].  Furthermore,  it  has 
been  shown  that  damage  growth  within  a  process  zone  during  fatigue  fracture  can  be 
described  by  a  self  -  similar  transformation  of  the  space  variables  [17].  Modeling  of 
such  fracture  processes  has  been  made  by  introducing  additional  kinematic 
parameters  on  the  basis  of  self  -  similarity  [16]. 

Formulations  of  constitutive  equations  for  damage  growth  and  criteria  for  crack 
initiation  are  well  recognized  by  now  to  be  important  problems  in  fracture  research. 
Although,  important  progress  has  been  made  in  this  direction,  a  number  of  questions 
still  remain  to  be  properly  addressed  and  resolved  for  a  better  understanding  of 
fracture.  In  this  two  Part  series  of  papers,  inquiries  into  the  characteristics  of  damage 
evolution  before  crack  initiation  and  criteria  for  crack  initiation  are  addressed.  In  Part  I, 
experimental  studies  on  damage  evolution  are  reported  using  amorphous  polystyrene 
as  a  model  material.  To  facilitate  experimental  observations,  the  site  of  crack  initiation 
is  located  in  space  by  inducing  a  60°  -V  notch  onto  the  mid-span  of  the  specimen 


2 


edge.  Hence,  efforts  are  concentrated  on  characterizing  damage  dissemination,  the 
associated  kinematic  parameters,  and  the  time  to  crack  initiation.  In  Part  II,  the  energy 
release  rate  due  to  damage  growth  is  evaluated  using  of  a  semi  -  empirical  method 
and  experimental  measurements.  On  the  basis  of  the  experimental  results  damage 
growth  in  front  of  the  notch  tip  is  modeled  with  a  first  order  reaction  equation. 

2.  EXPERIMENTAL  PROCEDURES 
MATERIAL  AND  SPECIMENS 

Commercially  available  amorphous  polystyrene  (PS)  from  Transilwrap,  Chicago,  IL, 
was  employed  in  the  studies  of  crack  initiation.  The  material  was  received  in  the  form 
of  sheets  with  dimensions  of  200x250  mm.  Strips  of  150x22x0.20  mm  were  cut  from 
these  sheets  with  a  razor  blade  and  sandwiched  between  two  pieces  of  acrylic  with 
dimensions  150x21x6.5  mm.  The  specimens  were  squared  up  to  150x20x0.20  mm 
with  the  use  of  a  fly  cutter  on  a  milling  machine.  The  speeds  of  the  fly  cutter  and  the 
longitudinal  power  feed  of  the  working  table  that  supports  the  specimens  were  11 
rotation -secT  and  0.23  mm  sec'',  respectively.  First,  sections  of  about  0.25  mm  thick 
were  cut  off  the  block  of  specimens  until  the  edges  of  the  specimens  were  on  the  same 
plane.  The  final  few  cuts  were  in  the  order  of  0.01  mm  so  that  the  induced  damage 
was  minimized.  The  block  was  then  flipped  over  and  the  same  procedure  was 
repeated  until  the  desired  width  was  achieved.  The  machined  edges  were 
metallographically  polished  to  a  0.5M.m  finish  to  prevent  formation  of  edge  crazes.  A 
double  angle  cutter  of  60°  was  used  to  notch  the  specimens  prior  to  their  removal  from 
the  block.  Subsequently  the  specimens  were  removed  from  the  block  and  washed 
carefully  with  distilled  water.  Note  that  this  procedure  of  specimen  preparation  ensures 
identical  notch  tip  geometries.  Finally,  the  specimens  were  annealed  in  a  temperature 
of  10°C  degrees  lower  than  Tg  (glass  transition  temperature)  for  48  hours  and  then 

allowed  to  slowly  cool  down  to  room  temperature.  Annealing  was  aimed  at  relieving 
any  residual  stresses  and  healing  any  damage  formed  during  the  cutting  procedures. 

EXPERIMENTAL  METHODS 

Tension  -  tension  fatigue  experiments  were  conducted  on  an  Instron  Testing  System 
in  laboratory  environment  at  ambient  temperature.  All  experiments  were  performed 
under  load  controlled  mode  with  sinusoidal  wave  form.  The  evolution  of  damage 
around  the  notch  tip  was  observed  by  means  of  a  traveling  optical  microscope 
attached  to  the  Instron  Testing  System.  The  fracture  process  was  recorded  using  a 
motor  driven  camera  that  was  attached  to  the  microscope.  Craze  distribution  was 
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evaluated  from  optical  micrographs  of  sectioned  specimens  (approximately  10-20|j.m 
thick)  which  were  prepared  by  standard  metallographic  and  polishing  procedures. 
Damage  density  distributions  were  obtained  by  covering  the  micrographs  of  the 
polished  sections  with  a  mesh  of  rectangles.  The  size  of  a  typical  mesh  depended 
upon  the  extent  of  damage  around  the  notch  tip  and  was  about  NixN2“  10x15  (Figure 
1).  In  all  cases  the  dimensions  of  the  rectangles  were  approximately  12x1 2pm.  In 
each  rectangle,  the  number  of  crazes  was  counted.  Craze  density  was  evaluated  as 
number  per  unit  area  p  [#/mm^]  or  as  p  [mm^/mm^].  p  represents  the  amount  of 

area  of  craze  mid  planes  per  unit  volume,  n  is  the  number  of  intersections  of  crazes 
with  the  vertical  test  line  at  the  respective  rectangle,  a  and  b  were  the  height  and  width 
of  a  rectangle,  respectively,  and  t  is  the  specimen  thickness. 

The  most  direct  way  to  experimentally  measure  damage  growth  rates  are  by 
examination  of  sectioned  samples  that  have  been  exposed  to  a  well  controlled 
loading  history.  In  this  way  damage  growth  is  "frozen  in"  at  different  degrees  of 
development.  Thus,  to  obtain  damage  distribution  and  growth  rates  before  initiation, 
the  experiments  designed  for  damage  evolution  studies  were  interrupted  at 
appropriate  time  intervals.  Subsequently,  the  specimens  were  polished  and  the 
number  of  crazes  in  each  rectangle  was  evaluated. 

The  approach  of  characterizing  damage  described  in  the  preceding  paragraphs 
implies  that,  at  a  given  time,  damage  distribution  around  the  stress  concentrator  is 
reproducible  for  a  certain  set  of  loading  conditions  and  specimen  geometry.  Damage 
accumulation,  however,  is  a  stochastic  process.  Accordingly,  damage  distribution 
should  result  from  averages  of  a  sufficiently  large  number  of  distributions  observed  on 
identical  specimens.  However,  such  an  effort  is  experimentally  difficult  to  carry  out. 
Commercially  available  PS  is  a  relatively  'ductile'  material.  This  results  in  a  large 
extent  of  crazing  during  crack  initiation  and  slow  crack  growth  [17,18].  Thus,  damage 
within  a  zone  in  this  material  is  amenable  to  a  statistical  treatment.  To  examine  the 
reproducibility  of  damage  distribution,  two  experiments  were  performed  under  the 
same  loading  conditions  and  specimen  geometry. 

The  time  to  crack  initiation  was  examined  under  the  loading  histories  shown  in 
Table  I.  Crack  initiation  was  considered  when  a  crack,  however  small,  would  appear 
at  the  notch  tip.  This  event  was  observed  with  the  use  of  a  Nikon  optical  microscope 
under  a  magnification  of  x200.  It  has  been  shown  that  crack  initiation  time  in  PS  under 
fatigue  loads  exhibits  a  variability  of  about  10%  [19].  This  scatter  was  considered 
small  and  thus,  efforts  were  concentrated  at  investigating  damage  growth  before 
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initiation  under  three  loading  conditions,  which  are  referred  to  as  C1 ,  C2,  and  C3. 
Four  experiments  were  conducted  under  Cl  and  three  under  each  of  C2  and  C3 
conditions.  The  loads  and  frequency  for  these  three  sets  of  experiments  are  shown  in 
Table  I. 

Due  to  difficulties  associated  with  measuring  an  accurate  crack  size  at  initiation 
using  an  optical  technique,  the  tests  were  interrupted  shortly  after  crack  initiation. 
Subsequently,  the  specimens  were  pulled  to  fracture.  A  crack  size  at  initiation  was 
assessed  by  examination  of  the  fracture  surfaces  and  its  dependence  on  the  stress 
level  was  evaluated  by  performing  fatigue  tests  under  different  loads. 

3.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 

Upon  application  of  the  load,  a  core  zone  of  crazes  appeared  at  the  notch  tip.  While 
this  zone  increased  in  size,  few  arched  crazes  nucleated  and  propagated  around  the 
core  zone.  As  time  progressed,  the  crazes  grew  longer  as  well  as  new  crazes 
appeared  in  and  around  the  vicinity  of  the  core.  A  series  of  optical  micrographs  taken 
during  a  typical  experiment  is  shown  in  Figure  2.  Crack  initiation  was  defined  at  the 
instance  when  a  crack  appeared  at  the  notch  tip. 

One  way  to  investigate  damage  evolution  is  to  compare  the  contours  of  equal 
damage  density  at  consecutive  configurations.  Such  pointwise  comparison,  however, 
may  not  be  appropriate  because  of  local  fluctuations  in  damage  density.  Indeed,  the 
contours  of  equal  damage  level  at  each  configuration  obtained  from  different 
specimens  exhibit  noticeable  fluctuation  on  the  scale  of  several  micrometers  (Figures 
3  to  7).  Therefore  integral  parameters,  such  as  the  moments  of  damage  distribution, 
were  used  here  to  evaluate  the  properties  of  damage  evolution  before  crack  initiation. 

Characterization  of  any  distribution  with  the  use  of  moments  should  involve 
comparison  of  a  sufficiently  large  number  of  moments.  However,  relatively  high  order 
moments  may  not  be  appropriate  to  compare  because  the  experimental  error  built  into 
the  measurement  makes  them  relatively  inaccurate.  In  the  present  analysis,  the 
distributions  of  damage  were  characterized  by  the  total  number  of  crazes  within  the 
zone  lo,  the  second  l2x.  Izy  and  fourth  Ux,  l4y,  and  l2x2y  central  moments.  These 
quantities  were  evaluated  from  the  following  expressions, 

Ni  N2  Ni  N2 

Jo  =  X  2  Pij  -  Imxny  =  X  £  [Xij-Xc]"'[yij-yc]"PijAyjAxi 

i=1  j=1  i=1  j=1 
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Here  i  =1  ,...Ni  and  j  =1  ,...N2  refer  to  the  size  of  the  mesh  (see  Figure  1 )  and  m,  n  stand 
for  the  order  of  the  moments  and  take  values  of,  0,  2,  or  4.  Xc  and  yc  are  the 
coordinates  of  the  damage  zone  center  with  reference  to  the  notch  tip  and  given  by, 

Ni  N2  Ni  N2 

Xc  =  X  Z  XiiPii  /  Z  £  Pij  -  Vc  =  0 

i=1  j=1  i=1  j=1 

yc  is  zero  due  to  symmetry  of  the  damage  zone  around  the  y  axis  (Figure  3). 

The  contours  of  equal  damage  density  p,  pertaining  to  two  specimens  fatigued 
under  conditions  C2  (Table  I)  and  up  to  the  instance  of  crack  initiation  are  displayed  in 
Figure  3.  Note  that,  the  contours  of  damage  density  have  been  symmetrized  with 
respect  to  the  bisector  of  the  V-notch  because  small  variations  in  density  were 
observed  at  some  points  located  at  equal  distances  from  the  X  axis.  These  differences 
were  attributed  to  local  fluctuation  in  damage  density.  Thus,  the  density  at  these 
locations  was  substituted  with  the  arithmetic  mean  of  the  respective  density  values. 
For  the  sake  of  visual  clarity,  the  level  of  experimental  error  is  not  shown  in  Figures  3 
to  6.  The  error  in  the  measurements  of  craze  density  far  from  the  notch  tip  was  about 
10%.  For  the  measurements  near  the  notch  tip,  the  error  was  about  15  -  20%.  The 
larger  value  of  error  was  because  the  high  craze  density  in  the  close  vicinity  of  the 
notch  tip  limits  optical  microscopy  which  makes  accurate  measurements  difficult.  From 
these  data  moments  of  the  distribution  with  respect  to  the  center  of  the  zone  as  well  as 
the  notch  tip  were  evaluated  (Table  II).  A  comparison  of  the  moments  shown  in  Table 
II  shows  that  the  largest  difference  is  ~10%  and  occurred  between  the  Ux  moments 
with  reference  to  the  notch  tip.  The  rest  of  the  differences  were  always  less  than  10%. 
These  data  and  the  reproducibility  of  damage  distribution  within  a  process  zone  in  the 
same  material  [17]  were  taken  as  sufficient  indications  that  damage  distributions 
during  crack  initiation  were  reproducible. 

Measurements  of  damage  density  within  the  zones  yielded  the  contours  of  equal 
density,  p,  that  are  shown  in  Figures  4,  5,  and  6,  for  conditions  Cl,  C2  and  C3, 
respectively.  From  the  experimental  data  presented  in  these  Figures  the  center  of  the 
damage  zone,  Xc,  can  be  calculated  at  each  configuration.  The  data  in  Figure  7 
indicate  that  the  center  of  gravity  of  the  damage  zone  increased  monotonically  with  the 
number  of  cycles.  To  investigate  the  type  of  transformation  the  damage  density  is 
undergoing  during  its  evolution,  the  following  characteristic  scales  along  the  X  and  Y 
axes  are  defined, 
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Since  I2X,  l2y  second  moments  of  the  damage  zone  and  Iq  is  the  zeroth 

moment  (which  expresses  the  total  amount  of  crazes  within  the  damage  zone)  and 
Oy  can  be  looked  upon  as  the  measurements  of  damage  spread  in  the  X  and  Y  axes, 
respectively.  Accordingly,  for  the  purpose  of  comparing  damage  distributions  within 
the  damage  zone  the  following  ratios  were  defined  between  ith  and  jth  configurations 
(i=  1, 2,  3  and  j=i+1). 
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h 
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The  type  of  transformation  of  damage  distribution  between  consecutive 
configurations  was  examined  by  comparing  the  following  quantities: 
(Xji,  Aji)),(iiji,  Mji)),  and  Hji)).  The  above  defined  ratios  for  conditions  C1,  C2, 
and  C3  are  shown  in  Table  III.  The  data  in  Table  III  indicate  a  relatively  large 
difference  in  the  ratios  between  the  first  and  second  configurations  for  conditions  C1. 
This  may  be  because  the  system  needs  some  time  for  damage  to  be  developed  so 
that  statistical  measurements  are  meaningful.  However,  differences  between  the  rest 
of  the  ratios  were  small. 

The  constancy  of  the  ratios  p-ji  between  configurations  shown  in  Table  III 

implied  that  the  evolution  of  damage  can  be  approximated  by  a  linear  transformation 
of  the  space  variables.  To  obtain  the  kinematic  parameters,  the  velocity  V^(x)  of  a 

point  within  the  zone  is  expressed  by  the  first  two  terms  of  Taylor  series  around  the 
center  of  the  damage  zone, 

Vm(x)  s=  Vnn(O)  +  Vnri,n(x)Xn  (3.1 ) 


Here,  the  first  term  represents  the  rigid  translation  of  the  damage  zone  and  the  second 
term  represents  the  rotation  (anti-symmetric  part  of  Vn,  „)  and  the  deformation 
(symmetric  part  of  „)  which  in  turn  can  be  decomposed  into  isotropic  expansion  and 
homogeneous  distortion  of  the  damage  zone.  In  this  section  the  velocity  Vn,(x)  is 
expressed  in  terms  of  the  transformation  coefficients  A.jj  and  Pjj.  Subsequently  the  rates 

of  expansion  and  distortion  are  evaluated. 

The  transformation  matrix  of  points  with  equal  damage  density  between 
configuration  i  and  j  is. 
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A,jj  0 


which  can  be  rewritten  as, 


1  +AX.jj  0 


0  1  +A|a.ji 


=  I+AD 


where  I  is  the  identity  matrix  and, 


AD  = 


AA,ji  0 
0  Ajjji 


Thus,  the  displacement  field  within  the  damage  zone  can  be  expressed  as, 

x+Axl  r  1+AA,ji  0  Irxi  r  AXc 

=  +  (3.5) 

_y+Ay  J  L  0  J  L  y  ■'  L  0  J 

where,  AXg,  represents  the  growth  of  the  center  along  the  x  axis.  Equation  (3.5)  can  be 
rewritten  as, 

r  Ax  1  [  A^-ii  0  Irxi  r  AXcl 


[  Ay  j  L  0  A)ijj  J  L  y  J  L  0  J 

Dividing  both  sides  by  the  increment  of  the  cycle  number  AN,  and  setting  >Vi, 

AN 

— >V2,  it  is  obtained, 

AN  AN  ^  AN  AN 

o][XKr>;cl 


V2J  [0  lijjiyi  10 


A  comparison  of  Equations  (3.1)  and  (3.7)  results  in. 


Vnn(x)=  ,  Vm(0)=[^  ,  Vm,n(x)  = 


Xji  0 
0  Ri. 
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Thus,  the  rates  of  expansion  e,  and  distortion  d,  of  the  damage  zone  can  be  expressed 
as, 


e  —  ^A,jj+}ijj),  and  d  —  (3.9) 

The  average  values  of  e  and  d  as  well  as  the  growth  of  the  gravity  center  of  the 
damage  zone  for  three  loading  conditions  are  displayed  in  Table  IV.  The  data  in  Table 
IV  show  that  while  the  growth  behavior  of  AX^j/AN  depended  upon  the  loading  history, 
the  average  rates  of  the  damage  zone  deformation,  Ae/AN  and  Ad/AN,  decreased 
until  crack  initiation. 

It  is  worth  pointing  out  that  damage  evolution  within  a  process  zone  on  the  same 
material  can  be  described  by  a  linear  transformation  [17].  These  findings  and  the 
results  of  the  present  studies  indicate  that  a  similar  transformation  relates  damage 
evolution  before  crack  initiation  and  during  slow  crack  growth. 

The  evolution  of  an  average  damage  density  <p>  =  E/A,  where  E  is  the  total 
number  of  crazes  within  the  zone  and  A  is  the  area  of  the  zone  with  normalized  cycle 
number  is  shown  in  Figure  8.  (Note  that  because  of  the  large  difference  in  the  cycles 
to  crack  initiation  between  the  different  experimental  conditions.  Cl;  300,000  cycles, 
C2:  30,000,  C3:  4,000,  the  number  of  cycles  has  been  normalized  with  the  respective 
cycle  number  at  initiation  in  order  to  clearly  see  the  trend  of  <p>).  It  is  interesting  to 
note  that  all  three  sets  of  data  can  be  approximated  by  an  equation  of  the  form  (lines  in 
Figure  8), 

<p>  =  1  -  exp(- KTN)  (3.10) 

where  N  is  the  number  of  cycles  and  X*  is  a  load  dependent  parameter  and  equal  to 
K  =  2.5x1 0‘5  (Cl),  8x10-5  (C2),  and  1.5x10-3  cycle-''(C3).  Expression  (3.10)  indicates 
that  the  dependency  of  <p>  on  the  cycle  number  is  of  a  damping  character.  Similar 
trends  of  damage  accumulation  have  been  reported  on  smooth  specimens  of  different 
polymers  under  creep  loads  [20]. 

The  data  in  Figure  8  show  that  while  the  rates  of  damage  accumulation  were 
different,  craze  densities  <p>,  at  crack  initiation  were  40,000,  44,200  and  47,000  mm’ 
respectively.  Note  that  the  error  between  the  extreme  values  of  these  densities  was 
16%  while  the  largest  difference  from  the  mean  was  about  10%.  Thus  one  may  be  led 
to  consider  the  value  of  <p>  as  being  independent  of  the  loading  conditions  at  crack 
initiation.  This  would  be  realistic  if  the  density  within  the  zone  was  constant  throughout 
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the  zone.  However,  the  results  of  these  studies  clearly  indicate  that  damage  was  not 
homogeneously  distributed  (Figures  3  to  6).  Thus,  an  average  quantity  like  <p>  may 
not  be  indicative  of  the  crack  initiation  event.  On  the  other  hand,  crack  initiation 
occurred  within  a  core  zone  ahead  of  the  notch  tip.  Morphologies  of  this  highly 
localized  damage  zone  at  initiation  for  conditions  C1  and  C3  are  shown  in  Figures  9 
and  10,  respectively.  Within  the  resolution  of  the  experimental  measurements, 
damage  density  within  the  core  was  about  1,300  mm^/mm^  and  approximately  the 
same  for  all  loading  conditions  (Figures  4  to  6).  Similar  value  of  damage  density  has 
been  obtained  in  the  close  proximity  of  the  crack  tip  grown  under  low  frequency  fatigue 
loads  [17].  The  constancy  suggests  that  a  certain  level  of  damage  density  is  required 
for  crack  initiation  and  subsequent  crack  growth  and  that  this  density  may  be  a 
material  parameter.  Additional  experimental  as  well  as  theoretical  research  is  needed 
to  firmly  establish  this  important  observation. 

The  optical  micrographs  in  Figures  9  and  10  bring  out  an  interesting  observation. 
While  the  density  at  the  notch  tip  was  practically  the  same,  the  cycle  numbers  to  crack 
initiation  were  drastically  different,  i.  e.,  300,000  and  4,000  cycles,  respectively. 
Moreover,  damage  dissemination  around  the  core  zone  was  quite  different  (Figures  9 
and  10).  Thus,  one  would  expect  that  in  the  case  where  large  crazes  surrounded  the 
core  zone  would  take  longer  to  initiate  a  crack  (Figure  10)  due  to  a  greater  extent  of 
shielding  of  the  peripheral  crazes  on  the  core.  However,  this  was  not  the  case. 
Therefore,  it  is  assumed  that  time  dependent  processes  within  the  damage  zone  play 
a  more  important  role  than  the  interaction  between  the  crazes. 

It  was  stated  earlier  that  a  small  crack  was  observed  at  the  notch  tip.  This  event 
was  associated  with  crack  initiation.  An  accurate  initial  crack  size,  however,  was 
difficult  to  detect  with  the  use  of  standard  optical  observations  during  the  fatigue  tests. 
To  evaluate  an  initial  crack  size,  the  specimens  were  pulled  to  fracture  shortly  after 
crack  initiation  and  the  fracture  surfaces  were  observed  under  an  optical  microscope. 
A  typical  micrograph  of  the  fracture  surface  morphology  near  the  notch  tip  is  displayed 
in  Figure  11.  Note  that  two  distinct  morphologies  can  be  seen  in  the  fracture  surface: 

(i)  a  relatively  rough  surface  from  the  notch  tip  which  extents  up  to  an  arched  front,  and 

(ii)  the  surface  with  tearing  like  features.  It  is  assumed  that  the  first  morphology  is 
associated  with  crack  initiation.  The  second  morphology  corresponds  to  the  crack 
growth  phase  resulting  from  the  simple  pulling  of  the  specimen  after  initiation.  The 
morphology  shown  in  Figure  11  suggested  that  the  crack  front  was  not  straight. 
Instead,  crack  grows  more  in  the  middle  than  at  the  edges  of  the  specimen  surface. 
The  crack  length  at  initiation  j( ,  taken  as  the  average  of  five  measurements  from  the 
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notch  tip  to  the  curved  front  is  shown  in  Figure  12  as  a  function  of  the  stress  ampiitude. 
These  data  indicate  that  the  dependence  of  the  crack  iength  on  the  applied  stress  is  of 
the  form,  Ln4  =  Bq  +  Bia  where  Bo  and  Bi  are  constants. 

Crack  initiation  times  Xj,  plotted  against  the  stress  amplitude  are  shown  in  Figure 
13  (note  that  Nj  =  VTj  where  Nj  is  the  number  of  cycles  for  crack  initiation  and  v  is  the 
test  frequency).  These  data  show  that  Xj  can  be  approximated  by  an  expression  of  the 
form,  Lnxj  =  Do  -  Dia  where  Do,  and  Di  are  parameters.  Similar  relationships  have 
been  obtained  for  different  materials  including  polymers  and  metals  under  constant 
loads  and  different  temperatures  [21 ,22]  and  have  been  explained  in  terms  of  stress  - 
temperature  activated  processes  at  the  molecular  level.  It  may  be  expected  that  such  a 
relation  would  not  apply  in  fatigue  loads.  However,  the  frequency  of  the  fatigue 
experiments  was  very  low.  Thus,  it  may  be  that  the  material's  response  is  primarily 
due  to  the  time  under  load  and  less  due  to  cycle  reversal.  Similar  results  were 
obtained  under  cyclic  loads  in  PANS  rocket  propellant  material  [23]  and  other 
polymeric  solids  [24]. 


TABLE  I:  Loading  Conditions  for  Crack  Initiation  Times 
and  Damage  Evolution  Studies 


^  max 

(MPa) 

min  (MPa) 

Frequency 

(sec-1) 

Cycles  to  Initiatio 
(x103) 

C1: 

8.00 

2.75 

0.33 

300±10 

C2: 

10.50 

2.75 

0.33 

30±2.5 

C3: 

14.30 

2.75 

0.33 

4±0.5 

12.00 

3.00 

0.33 

20±2.0 

13.00 

3.00 

0.33 

10±1.2 

TABLE  II:  Inertia  Moments  for  two  Configurations  of  Damage  before  Crack  Initiation 


PWWBBHW 

1 

2 

1 

2 

Reference 

Notch  Tip 

Notch  Tip 

Damage  Zone 

Damage  Zone 

Center 

Center 

Xc(^m) 

55.0 

52.0 

Yc 

0.0 

0.0 

Iq  (number) 

732 

680 

lac  (mm2) 

2.8 

2.5 

0.65 

0.63 

Igv  (mm2) 

0.46 

0.42 

0.46 

0.42 

Ux  (mnfi^) 

0.20x10-1 

0.18x10-1 

0.11x10-2 

0.11x10-2 

I4V  (mm^) 

0.57x10-3 

0.52x10-3 

0.57x10-3 

0.52x10-3 

I2x2v(mni^) 

0.21x10-2 

0.18x10-2 

0.40x10-3 

0.39x10-3 

12 


TABLE  III:  Moment  Ratio  of  Damage  Distribution  at  Consecutive  Configurations 
before  Crack  Initiation 


Condition  Cl 

^21=1. 20 

|l2i=1-02 

X21  P2i=1.22 

X32=1.57 

EBB 

|l32=1 .30 

IStSEBI 

^2  ll32~2.04 

X.43=1.43 

H43=1.30 

M43=1.30 

^43  1^43=1  -®6 

Hiq=1.90 

Condition  C2 

X21=1.11 

A2i=1.10 

1121=1-22 

A^1 1X21  -35 

BSISH 

>92=1 .17 

EBm 

|l32=1  -23 

^2  M32“1 

EBN 

Condition  C3 

X2i=1.60 

A2i=1.58 

p.2i=1.38 

dSEEH 

M^21=2-21 

^32=1 .45 

1132=1-36 

^2  M32=1  -37 

TABLE  IV:  Growth  Rates  of  the  Gravity  Center,  Isotropic  Expansion  and  Distortion  for 
three  Loading  Conditions 


Configuration 

AXc/AN  (x104) 
(pm/cycle) 

Ae/AN  (x106) 
1/cycle 

Ad/AN  (x106) 

1 /cycle 

Conditions  C1 : 

1-2 

1.3 

13.0 

4.4 

2-3 

1.5 

3.0 

1.0 

3-4 

1.6 

2.0 

0.7 

Conditions  C2: 

1-2 

4.5 

20.0 

4.5 

2-3 

11.8 

15.0 

3.6 

Conditions  C3: 

1-2 

141.5 

27.0 

48.0 

2-3 

4.5 

1.75 

3.0 
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Fig.  3  Contours  of  equal  damage  density  p  [mm2/mm3],  for  two 

specimens  fatigued  under  the  same  conditions  and  equal  times 
(Conditions  C2) 
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Fig.  4  Contours  of  equal  damage  density  p  [mm^/mm^],  at  four 

consecutive  configurations  before  crack  initiation  (Conditions  C1) 
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Fig.  5  Contours  of  equal  damage  density  p  [mm2/mm3],  at  three 

consecutive  configurations  before  crack  initiation  (Conditions  C2) 
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Contours  of  equal  damage  density  p  [mm2/mm3],  at  three 
consecutive  configurations  before  crack  initiation  (Conditions  C3) 


Fig.  7  Evolution  of  the  gravity  center  of  the  damage  zone 
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Fig.  8  Normalized  average  density  before  crack  initiation  plotted  against 
normalizes  cycle  number 
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Fig.  9  Optical  micrographs  showing  damage  distribution  at  crack 
initiation  (Conditions  C1). 
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10  Optical  micrographs  showing  damage  distribution  at  crack 
initiation  (Conditions  C3) 
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4  ANALYSIS 

ENERGY  EVALUATIONS 

The  experimental  results  reported  in  the  last  section  [25]  demonstrated  that  during  the 
process  of  crack  initiation  in  amorphous  polystyrene  (PS)  two  distinct  craze  patterns 
appeared  around  the  stress  concentrator;  a  core  zone  of  dense  crazes  iocated  in  front 
of  the  notch  tip  and  peripheral,  less  dense  crazes,  surrounding  the  core.  A  typical 
schematic  of  damage  before  crack  initiation  is  shown  in  Figure  14.  While  the  extent  of 
the  peripheral  crazes  depended  on  the  loading  history,  the  density  within  the  core 
zone  at  initiation  was  found,  to  a  good  degree  of  approximation,  independent  of  the 
loading  conditions.  Moreover,  it  was  shown  that  damage  evolution  before  crack 
initiation  can  be  described  by  a  self  -  similar  transformation.  These  results  allowed  to 
deduce  the  movements  of  the  damage  zone  to,  translation,  isotropic  expansion,  and 
homogeneous  distortion  [25].  Although,  these  results  provided  important  information 
in  the  efforts  to  characterize  the  crack  initiation  process,  a  further  understanding 
requires  calculation  of  the  energy  release  rates  associated  with  these  kinematic 
parameters.  Knowledge  of  the  energy  release  rates  would  allow  to  formulate  kinetic 
equations  for  damage  growth. 

Calculation  of  the  stress  and  strain  fields  in  an  elastic  body  with  defects  is  the  first 
step  towards  the  evaluation  of  the  energy  release  rates  or  stress  intensity  factors. 
However,  such  a  task  is  not  always  easy  to  accomplish.  The  complex  interaction 
between  the  damage  elements  and  the  crack  makes  an  analytical  solution  very  difficult 
to  obtain.  Therefore,  various  approximations  have  been  proposed  in  the  literature.  In 
works  reported  in  [26-28],  damage  is  modeled  by  an  effective  inclusion.  The 
properties  of  the  inclusion,  however,  are  not  easily  determined  experimentally.  In 
addition,  a  substitution  of  an  array  of  defects  by  a  homogeneous  medium  may  not  be 
realistic  when  defect  distribution  within  a  damage  zone  is  heterogeneous.  In 
references  [29-33],  solutions  for  the  stress  and  strain  fields  are  based  on  a  detailed 
description  of  the  defect  size,  location,  and  orientation.  The  mathematical  difficulties 
involved  in  such  a  description  are  enormous  and  solutions  to  certain  configurations  of 
crack  -  microcrack  systems  are  possible.  In  an  approach  reported  in  [34],  a  zone  of 
defects  around  a  crack  is  characterized  by  statistical  distributions  of  defect  density, 
size,  and  orientation.  Average  parameters  are  then  introduced  to  model  the  effects  of 
a  damage  zone  on  energy  release  rates  or  stress  intensity  factors.  However,  when  a 
large  number  of  defects  is  present,  it  is  difficult  to  produce  experimentally  accurate 
distributions.  A  semi  -  empirical  analysis  has  also  been  proposed  to  evaluate  an 
energy  release  rate  due  to  damage  growth  [35].  It  is  based  on  crack  or  craze  openings 
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and  the  double  layer  potential  technique.  This  method  was  used  in  the  present  work 
as  a  first  approximation  to  estimate  energy  changes  due  to  damage  using  the 
experimental  data  and  certain  assumptions  regarding  craze  openings. 

An  important  aspect  of  the  initiation  process  is  the  condition  under  which  crack 
initiation  takes  place.  Although  various  mechanisms  have  been  proposed  in  the 
literature  to  explain  crack  initiation  in  different  materials,  quantitative  descriptions  of 
initiation  have  been  difficult  to  formulate.  In  a  physically  sound  model  presented  in 
[36]  crack  initiation  is  taken  as  the  instance  when  the  sum  of  the  displacements  of  the 
accumulated  dipoles  before  initiation  exceeds  the  theoretical  displacement  to  cause 
fracture  in  a  perfect  crystal.  The  model  derives  expressions  for  the  time  to  crack 
initiation  and  the  length  of  the  initiated  crack  and  is  in  good  agreement  with 
experimental  results.  It  is  not  clear,  however,  if  the  model  would  apply  to  different 
modes  of  fatigue  damage.  Principles  of  continuum  damage  mechanics  have  also 
been  employed  for  crack  initiation.  In  these  models,  it  is  assumed  that  a  macrocrack 
appears  when  the  elastic  stiffness  drops  to  zero,  or  when  a  damage  parameter 
reaches  a  critical  value  [37,  38].  However,  damage  parameters  employed  in 
continuum  damage  mechanics  may  not  always  be  realistic.  Models  based  on 
irreversible  thermodynamics  have  also  been  reported,  aithough,  they  have  not  been 
tested  experimentally  [39,  40].  Constitutive  models  for  nucleation,  growth,  and 
coalescence  of  various  submicroscopic  failure  processes  have  also  been  proposed 
(for  a  review  see  [41]).  These  models  are  of  two  types,  (a)  failure  initiation  is  activated 
by  stress  and/or  temperature,  and  (b)  failure  initiation  is  due  to  plastic  deformation.  In 
the  stress  -  temperature  models,  the  effects  of  a  tensile  stress  on  the  free  energy  for 
thermally  activated  microscopic  processes,  such  as,  atomic  bond  breaking  and/or 
diffusion  of  vacancies  or  polymer  chain  segments  are  considered.  The  plastic 
deformation  models  are  primarily  based  on  heterogeneous  slip  due  to  a  shear  stress. 
Typical  examples  are  grain  boundary  sliding,  and  plastic  deformation  around 
inclusions.  Thermally  activated  processes  may  also  play  a  role  in  the  plastic 
deformation  models.  This  is  when  the  initially  static  dislocation  or  other  defects 
become  mobile  upon  application  of  stresses  and  the  associated  changes  in 
temperature.  These  models  describe  well  certain  aspects  of  submicroscopic  flaw 
nucleation  and  have  been  used  in  computer  simulations  with  varying  degrees  of 
success  [41 , 42]. 

In  this  paper,  a  semi  -  empirical  method  and  the  experimental  results  reported  in 
Sec.  3  are  used  to  evaluate  an  energy  release  rate  associated  with  damage  growth 
before  crack  initiation.  Subsequently,  the  rates  of  the  damage  zone  movements  are 
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correlated  with  the  energy  release  rate.  Next,  attention  is  focused  on  modeling  the 
evolution  of  damage  within  the  core  zone  (Figure  14).  Guided  by  the  experimental 
results  [25]  and  works  by  various  researchers  on  fracture  of  polymeric  solids  [22,  43- 
51],  it  is  assumed  that  the  particular  fracture  process  reported  in  [25]  is  due  to 
submicroscopic  flaw  nucleation  resulting  from  thermal  fluctuations.  Thus,  the  ideas  of 
fracture  kinetics  are  employed  to  characterize  damage  evolution  at  the  core  zone  with 
a  first  order  kinetic  equation. 

The  experimental  results  presented  in  Part  1  of  these  investigations  indicated  that 
damage  evolution  before  crack  initiation  can  be  described  by  translation  and 
deformation  of  the  damage  zone.  Thus,  the  changes  in  the  total  potential  energy  of  the 
system  can  be  expressed  as. 


-  APt  =  Ji AXc  +  MAe  -I-  NjjAdij 

Here,  AXc,  Ae  and  Adjj  are  the  increments  of  the  damage  zone  center,  isotropic 
expansion,  and  distortion,  respectively,  and  Ji,  M,  and  Njj  are  the  corresponding 
enrage  release  rates.  To  investigate  the  evolution  of  the  damage  zone  before  crack 
initiation  these  energy  release  rates  should  be  evaluated.  Efforts  to  calculate  these 
energies  have  been  reported  in  the  literature  [34].  However,  use  of  these  methods 
requires  an  accurate  description  of  the  craze  length,  their  distribution,  and  interaction. 
Such  complete  characterization  is  difficult  to  carry  out  in  systems  where  a  dense 
damage  zone  precedes  crack  initiation.  In  this  work  a  total  energy  release  rate  A, 
defined  as. 


(4.1) 


is  evaluated  using  a  semi  -  empirical  analysis  and  experimental  measurements.  Here 
4  is  a  characteristic  size  of  the  damage  zone,  /q  is  the  notch  depth  and  t  is  the 
specimen  thickness.  With  respect  to  the  characteristic  size,  various  length  parameters 
of  the  damage  distribution  can  be  used.  For  simplicity,  4  is  taken  as  the  horizontal 
distance  of  the  outermost  damage  contour  from  the  notch  tip  [1].  Such  a  quantity  is 
easily  evaluated  from  the  experiments  and  due  to  self  -  similarity  of  damage  it 
accounts  for  the  entire  size  of  the  damage  zone. 

The  total  potential  energy  of  solid  with  a  crack  (or  a  notch)  and  an  array  of 
microdefects  is, 

P|-Pc*P<i*Pi*Po  (4.2) 


29 


The  first  two  terms  of  the  right  hand  side  are  the  potential  energies  associated  with  the 
main  crack  and  the  array  of  microdefects,  respectively,  P|  is  the  potential  energy  due  to 
interaction  of  the  crack  and  microdefects,  and  Pothe  potential  energy  of  the  initial  state 
(i.  e.,  without  a  crack  and  mirodefects).  Each  term  in  (4.2)  can  be  defined  as  the 
difference  between  the  strain  energy  of  the  solid  and  the  work  done  by  the  applied 
load.  For  a  linear  elastic  solid,  the  total  potential  energy  can  be  written  as, 

p,=  -If  A,  (4.3) 

where  F  is  the  boundary  load  which  is  taken  as,  F  =  amaxA  (ctmax 's  the  maximum 
stress  of  the  fatigue  cycle,  and  A  is  the  specimen's  cross  section)  and  At  is  the  total 
displacement  at  the  points  of  load  application.  Since  the  load  is  prescribed,  the  total 
displacement  at  the  specimen  boundaries  are  evaluated  using  a  semi  -  empirical 
method  for  energy  calculations  [35].  To  carry  out  this  task,  the  following  simplifications 
and  assumptions  are  employed, 

1 .  Plane  stress  elastic  deformations  are  assumed  throughout  the  specimen  and 
within  the  damage  zone. 

2.  Residual  stresses  are  negligible. 

3.  The  principle  of  superposition  applies. 

4.  The  Green's  tensor  for  the  displacement  a  point  X  due  to  a  unit  discontinuity  at  ^  is 
not  available  for  a  semi  -  infinite  plate.  Thus,  the  displacements  at  the  grips  of  the 
specimen  are  obtained  by  using  the  Green's  tensor  for  an  infinite  plate.  This  is  only  an 
approximation  since  the  damage  zone  is  very  close  to  the  notch  boundaries  (Figure 
14)  and  thus,  the  stress  and  strain  fields  are  influenced  by  these  boundaries. 

To  estimate  the  extent  of  the  error  caused  by  the  specimen  boundary,  the 
displacements  at  the  grips  of  the  specimen  were  calculated  for  the  cases  when  a  crack 
equal  to  the  notch  depth  is  in  middle  of  the  specimen  and  for  a  specimen  with  an 
edged  crack  of  the  same  length  [52].  The  displacements  at  the  grips  of  the  single  edge 
cracked  specimen  were  found  greater  by  8%.  Thus,  in  the  present  analysis,  the 
displacements  calculated  using  the  Green's  tensor  for  an  infinite  plate  were  multiplied 
by  1.08. 

The  principle  of  superposition  applied  to  the  problem  at  hand  results  in  the  total 
displacement, 

At=Aeiastic+Acrack+Acrazes  (4"  4) 
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Next,  each  term  in  Equation  (4.4)  should  be  evaluated.  It  was  argued  earlier  that 
during  the  process  of  crack  initiation,  crack  growth  from  the  notch  does  not  occur. 
Moreover,  it  is  assumed  that  the  deformation  due  to  the  bulk  of  the  specimen  does  not 
change  between  configurations.  Accordingly,  the  only  contribution  to  the  energy 
release  rate  A,  is  that  due  to  damage  growth.  With  reference  to  Figure  15a,  the 
displacement  vector  Ui(x,^)  at  point  x  due  to  a  discontinuity  opening  bj(^)  is  given  by, 

Ui(X,5)  =  |  b,(^)4.ij(x4)d4  (4.5) 

Here  <^jj(x,4)  is  the  second  Green’s  tensor  which  stands  for  the  displacement  vector  at 
point  X  due  to  the  unit  discontinuity  at  4. 

(niR,-n|R,+n,Rkdii)  +  ^  RkRiR,]  (4.6) 

where  v  is  the  Poisson's  ration,  is  the  unit  normal  to  the  discontinuity  and  Ri=  Xj-^ 

andR2==R?  +  R|.  The  integration  is  performed  over  the  length  of  the  discontinuity,  s. 
Employing  the  coordinate  system  shown  in  Figure  15,  and  assuming  that  on  the 
average  sense  crazes  are  horizontal  at  any  point  within  every  square  of  the  grid 
(Figure  15),  ni=0,  n2=1,  and  Ri=(xr^i),  R2=  (X2  -  ^2)- 

To  continue,  the  shape  of  the  discontinuity  b(^)  should  be  identified.  With  reference 
to  the  Figure  15,  the  crazes  within  a  square  are  taken  horizontal  and  segregated  to  a 
single  discontinuity  of  length  equal  to  the  length  of  the  square.  The  width  of  the 
discontinuity  is  assumed  to  be  equal  to  the  product  of  the  some  reference  opening  Pq. 
and  the  number  of  crazes  within  the  respective  square.  Thus,  at  the  typical  square  of 
ith  column  and  jth  row  with  a  number  of  crazes  Py  which  is  measured  experimentally 
the  opening  of  an  effective  discontinuity  is, 

b|f(^)=toP,)(y  (4.7) 

Accordingly,  the  vertical  displacement  at  a  point  X(xi,  X2)  due  to  the  effective 
discontinuity  within  a  square  whose  center  is  located  at  ^(|i,  ^2)  's  obtained  by 
combining  Equations  (4.5)  to  (4.7), 


where  c  is  the  size  of  the  square  in  the  grid  (Figure  15b).  The  dispiacement  at  the  grips 
of  the  specimen  due  to  a  discontinuity  was  obtained  by  a  numericai  integration  of  (4.8) 
using  a  standard  routine,  at  X2=h,  where  h  is  half  the  specimen  height. 

The  total  dispiacement  at  the  grips  is  obtained  by  summing  up  the  contributions  of 
the  individual  displacements  caused  by  each  of  the  effective  discontinuities  within  the 
array, 


Ni  N2 

U2(xi,h)  =  X£  U‘|(xi,h)  (4.9) 

1=1  i=i 

The  average  displacement  at  the  grips  due  to  the  damage  zone  Acraze.  is  obtained  by 
averaging  the  displacements  across  the  specimen  width  B, 

^raze  =  U2(Xi,h)dXi  (4.10) 

A  parameter  that  remains  to  be  identified  is  a  reference  craze  opening  bg  (Equation 
4.7).  Due  to  difficulties  in  evaluating  such  a  quantity  in  the  experimental  studies,  an 
opening  of  bo  =  was  assumed.  It  can  be  seen  from  Equation  (4.7)  that  the  results 
are  directly  proportional  to  bo .  Thus,  the  trends  of  the  energy  release  rate  remain 
unaltered. 

Values  of  the  characteristic  size  4,  cycle  number,  and  energy  release  rates  A, 
(Equation  4.1)  for  two  loading  conditions  Cl  and  C2  are  shown  in  Table  V.  Using 
these  energy  values,  the  evolution  of  the  damage  zone  size,  A4/AN,  isotropic 
expansion,  Ae/AN,  and  distortion,  Ad/AN  [1],  are  shown  in  Figures  16  (conditions  Cl) 
and  4  (conditions  C2).  The  data  in  these  Figures  suggest  that  the  overall  growth 
behavior  of  the  damage  zone  before  crack  initiation  decreases  monotonically.  This 
trend  is  consistent  with  the  growth  behavior  of  an  average  damage  density  before 
crack  initiation  (sec.  3).  Moreover,  it  is  interesting  that  the  monotonic  decrease  shown 
in  Figures  16  and  17  resembles  the  behavior  of  small  cracks  observed  in  different 
materials  [53]. 

The  correlations  shown  in  Figures  16  and  17,  although  approximate,  illustrate  the 
trends  of  damage  growth  before  crack  initiation.  However,  criteria  for  crack  initiation 
can  not  be  easily  formulated  in  terms  of  some  Griffith’s  like  criteria.  That  is,  for  the 
fracture  process  investigated  in  this  work,  an  energy  release  rate  (or  stress  intensity 
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factor)  quantity  associated  with  crack  initiation  would  depend  on  the  specific  energy  for 
craze  formation,  the  amount  of  crazes  as  a  function  of  time  and  applied  loads,  as  well 
as  the  interactions  between  the  crazes  and  the  stress  concentrator.  It  is  difficult, 
however,  to  evaluate  these  contributions  with  the  existing  experimental  and  analytical 
methods.  Answers  to  these  questions  should  await  the  development  of  more  refined 
analytical  methods  and  improved  experimental  techniques. 

The  experimental  results  of  this  work  [25]  have  shown  that  damage  before  crack 
initiation  consisted  of  a  core  of  dense  crazing  and  a  peripheral  less  dense  zone  of 
crazing.  While  the  pattern  of  the  peripheral  crazes  depended  upon  the  applied  load, 
the  density  within  the  core  zone  was  independent  of  the  loading  conditions.  In  the 
remaining  of  the  paper,  this  later  observation  is  combined  with  principles  based  on 
microstructural  models  for  defect  nucleation  and  certain  assumptions  to  describe 
damage  accumulation  in  the  core  zone  and  the  instance  of  crack  initiation. 

KINETIC  CHARACTERISTICS  OF  DAMAGE  IN  THE  CORE  ZONE 
Experimental  and  analytical  investigations  have  shown  that  fracture  of  polymeric 
materials  is  due  to  thermal  fluctuations  on  the  molecular  level.  Thus,  reaction  rate 
theory  has  been  employed  to  model  their  behavior.  The  exact  mechanisms  that  lead  to 
failure  of  these  materials,  however,  are  not  known  a  priori.  Consequently,  various 
hypotheses  have  been  proposed.  In  one  of  them,  it  is  assumed  that  the  primary 
molecular  mechanism  giving  rise  to  fracture  is  scission  of  polymer  molecules  under 
stress.  In  such  a  case,  an  activation  energy  is  related  to  the  energy  of  thermal 
destruction  [22,  45,  49].  Another  hypothesis  is  based  on  polymer  chain  slippage  and 
rearrangement  upon  stress  [51].  For  such  processes  an  energy  barrier  is  determined 
by  the  rate  of  polymer  molecules  slipping  over  each  other.  A  combination  of  bond 
breakage  and  slippage  has  also  been  proposed  [49,  51].  The  extent  to  which  either  of 
these  two  mechanisms  dominates  depends  primarily  on  material  and  test  temperature. 
However,  it  is  difficult  to  quantify  their  contributions.  Thus,  they  are  usually  treated  as  a 
mechanism  with  a  single  energy  barrier.  Whatever  the  exact  mechanisms  may  be,  it 
appears  that  models  based  on  reaction  rate  theory  describe  certain  features  of 
polymer  fracture  reasonably  well,  although,  certain  modifications  and  assumptions 
may  be  necessary  to  accommodate  experimental  observations  [48]. 

To  model  the  process  of  damage  accumulation  within  the  core  zone  (Figure  14),  it 
is  assumed  that  on  the  submicroscopic  level  fracture  is  the  result  of  thermal 
fluctuations  of  polymer  chain  segments  which  form  crazes.  Within  the  crazes  chain 
rupture  and/or  chain  pull  -  out  occur  which  reduce  the  load  bearing  capacity  of  the 
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fibrils  thus,  causing  additional  ruptures  or  chain  pull  -  out.  This  is  the  stage  of 
submicrocrack  development.  When  the  number  of  submicrocracks  reaches  a  critical 
level,  a  macroscopic  crack  appears. 

To  describe  this  process  in  quantitative  terms,  it  is  assumed  that  there  exists  a 
certain  number  of  elements  Rq  within  the  area  around  the  notch  tip  where 
submicrocracks  may  nucleate  [47].  Therefore,  the  rate  of  submicrocrak  growth 
dR(t)/dt,  is  postulated  according  to  the  absolute  rate  theory, 

^=Kb(Ro-R(t))  (4.11) 

where  [seo-'']  can  be  looked  upon  as  the  rate  of  submicrocrack  accumulation.  It 
should  be  noted  that  the  contribution  of  defect  healing  is  not  taken  into  account  in  the 
rate  equation.  It  is  shown  later  in  this  work  that  its  contribution  is  much  smaller  in 
comparison  to  that  due  to  submicrocrak  nucleation.  Thus,  it  can  be  safely  neglected  in 
the  present  analysis. 

In  the  experimental  studies,  the  extent  of  crazes  in  the  core  zone  was  measured  at 
different  times  before  crack  initiation.  Therefore,  it  is  assumed  that  the  amount  of 
submicrocracks  at  any  time  prior  to  crack  initiation  is  proportional  to  the  extent  of 
crazes  within  the  core  zone.  Furthermore,  having  in  mind  that  damage  density  within 
the  core  reached  a  constant  value  [25],  independent  of  the  loading  conditions,  a 


normalized  damage  parameter  is  introduced  in  the  following  form, 

r(t)=-^  (4.12) 

ric 

where  Rg  is  the  density  of  submicrocracks  at  initiation.  On  the  basis  of  this  definition 
r(t)  takes  the  following  limiting  values, 

r(t=0)  =  0  corresponds  to  no  damage  (4.1 3a) 

r(t=Xi)  =  1  corresponds  to  crack  formation  (4.1 3b) 

where  Tj  is  the  crack  initiation  time.  Accordingly,  Equation  (11)  may  be  written  as, 

^  =  -Mt)  +  Kbro  (4.14) 


where  the  ro=Ro/Rc  is  a  constant. 
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Assuming  that  the  submicroscopic  mechanisms  taking  piace  during  damage 
accumuiation  are  driven  by  thermai  fluctuations  and  the  applied  stress,  the  rate 
constants  for  nucieation  and  healing  are  given  by  [43], 

Kb=J-exp[-{Uo-AU/2))/kT]  (4.15a) 

Kh  =;^xp[-(Uo+AU/2)/kT]  (4.1 5b) 

T'O 

Here  1/xo  is  the  vibration  frequency,  which  for  most  solids  it  is  found  to  be  in  the  order 
of  10'*3sec’'*;  Uq  [J/mole]  is  the  activation  energy  due  to  thermal  fluctuations;  AU 
[J/mole]  is  the  change  in  activation  energy  due  to  the  applied  stress;  k=1 .38X1 0*21 
[J/K]  is  the  Boltzmann's  constant;  T  [K]  is  the  temperature.  The  terms 
exp[-(Uo-AU/2))/kT]  and  exp[-(Uo+AU/2))/kT]  represent  the  probabilities  that  a 
weak  region  is  at  an  energy  level  Uo-AU  or  Uo+AU. 

To  estimate  the  rate  constant  from  the  above  expressions  both  energy  quantities 
Uq  and  AU  should  be  evaluated.  One  way  to  achieve  this  goal  is  to  use  the  methods 
of  statistical  physics  that  would  require  knowing  enough  about  the  types  of 
microdefects,  their  distribution  in  space  and  time,  and  the  range  of  their  interaction. 
Such  an  effort,  however,  is  very  difficult  to  achieve  for  practical  systems.  To  circumvent 
this  problem,  the  dependence  of  the  energy  barrier  on  stress  is  usually  approximated 
with  potentials  appropriate  to  a  particular  molecular  mechanism  of  failure  [54].  In  this 
work  we  resort  to  an  empirical  method  for  evaluating  Uq  and  AU  that  has  been  shown 
to  apply  well  in  different  materials. 

Extensive  experimental  research  on  various  solids  subjected  to  static  or  cyclic 
tensile  stress  o,  have  demonstrated  that  the  time  to  fracture  Xj,  can  be  described  by 
[22,55], 

Ti  =  xbexp(y^)  (4.16) 

where  y  [J/mole-MPa]  is  a  parameter  that  describes  the  degree  of  non  uniformity  of 
stress  distribution  in  the  loaded  specimen.  The  magnitude  of  Uq  has  been  found  to 
coincide  with  the  interatomic  binding  energy  of  various  materials  [22].  Both  quantities 
Uq  and  y  are  obtained  by  fitting  Equation  (4.16)  to  experimental  results  of  time  to 
fracture  where  c  is  taken  as  the  applied  stress  in  creep  experiments  and  as  the  stress 
amplitude  in  the  case  of  fatigue  tests.  Evidently,  the  time  to  fracture  given  by  Equation 
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(4.16)  is  equal  to  the  inverse  of  the  rate  constant  in  Equation  (4.15a)  when  the 
reduction  in  the  activation  energy  is  linearly  related  to  the  stress,  i.  e.,  AU  =  2'yj. 

Although,  Equation  (4.16)  describes  failure  times  of  smooth  specimens  under 
constant  and  fatigue  loads  very  well,  it  has  not  been  applied  to  crack  initiation  times 
from  a  stress  concentrator.  This  may  be  because  in  smooth  specimens  the  applied 
stress  field  is  assumed  uniform  and  thus,  satisfies  the  requirements  of  the  equilibrium 
Boltzmann’s  statistics.  On  the  other  hand,  the  high  stress  gradients  around  a  stress 
concentrator,  presumably  violate  the  requirements  of  Boltzmann's  statistics. 
Accordingly,  it  may  be  physically  unrealistic  at  first  to  apply  (4.16)  to  data  of  crack 
initiation  times  from  a  stress  concentrator.  The  experimental  results  have  shown  a 
relatively  large  number  of  crazes  at  the  core  zone  (Figure  14),  and  that  this  density 
reaches  a  load  independent  level  at  the  instance  of  crack  initiation  (sec.  3).  Thus,  it  is 
assumed  that  the  stress  level  over  the  core  zone  (Figure  14)  immediately  ahead  of  the 
notch  tip  is  uniform  and  proportional  to  the  amplitude  to  the  applied  stress,  in  terms  of 
fracture  mechanics,  this  is  similar  to  considering  the  layer  of  damage  as  a  cohesive 
zone  similar  to  that  in  the  Dugdale  -  Barenblatt  model  [56]. 

To  examine  the  crack  initiation  times  according  to  (4.16)  the  stress  parameter  is 
approximated  by  the  level  of  mean  stress,  aa=(cymax-<Jmin)-  Under  the  assumptions 
stated  in  the  preceding  paragraphs,  parameters  Uq  and  y  are  evaluated  by  fitting  the 
experimental  data  on  crack  initiation  reported  in  [25].  For  this  purpose.  Equation  (4.16) 
is  written  as, 

CTa  =  -iiLnTi+J^LnTo+^  (4.17) 


Using  the  data  in  Figure  13  [1],  a  plot  of  Oa  against  LnXj  for  five  experiments  resulted 
in,  ^  =  1.46  and  ^  Ln  ^0+^  =  25.20.  Solving  for  y  and  Uq  one  obtains  y  =  2.08 


kJ/MPa  mole  and  Uq  =  1 17.13  kJ/mole. 

It  has  been  reported  that  for  PS,  the  energy  for  thermal  destruction  is  about  230 
kJ/mole  [22]  and  an  activation  enthalpy  for  secondary  chain  motions  below  glass 
transition  temperature  about  65  kJ/mole  [34].  If  the  activation  energy  indicates  the 
types  of  submicroscopic  mechanisms,  a  value  of  Uq  =  115.6  kJ/mole  suggests  that 
secondary  chain  motions  as  well  as  chain  scission  are  possible  rate  determining 
processes  which  lead  to  crack  initiation.  This  is  consistent  with  the  observation  that 
dense  crazing  precedes  crack  initiation. 
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The  preceding  analysis  of  the  experimental  data  on  crack  initiation  times  has 
resulted  in  the  evaluation  of  an  activation  energy  Uoand  the  changes  due  to  the 
applied  stress  AU  =  27^3-  Therefore,  the  rate  constants  Kb  and  K^,  can  be  calculated. 
Using  Equations  (15)  and  the  values  of  Uq  and  y  obtained  from  the  foregoing  analysis 
it  is  obtained,  Kt/K^  ~  10^.  That  is,  healing  processes  may  be  neglected  in  the  rate  of 
submicrocrack  development.  Accordingly,  since  Kb  is  independent  of  time,  the 
solution  of  Equation  (4.14)  is, 

r(t)  =  ro  +  Cexp(-Kbt)  (4.1 8) 

Imposing  the  condition,  r(0)  =  0,  one  obtains, 

r(t)  =  ro[1  -exp(-Kbt)]  (4.19) 

The  value  of  Tq  is  estimated  from  the  condition  r(t=Xj)  =  1 .  For  conditions  Cl  an  C2,  it 
was  fond  that  Pq  =  1 .58,  while  for  C2,  Pq  =  2.00.  This  larger  value  was  attributed  to  the 
fact  that  the  crack  initiation  time  under  conditions  C2  was  about  40%  smaller  than  that 
predicted  by  (4.16).  The  solid  lines  in  Figure  18  represent  Equation  (4.19)  and  the 
points  are  the  experimental  data  of  the  normalized  density  for  loading  conditions 
Cl  (a),  C2(b)  and  C3(c),  respectively.  Because  of  the  large  differences  in  the  crack 
initiation  times,  the  data  for  each  condition  are  plotted  separately  for  the  sake  of  visual 
clarity.  The  data  in  Figure  18  demonstrate  that  Equation  (4.19)  fits  the  data  reasonably 
well.  Expressions  similar  to  (4.19)  have  been  also  reported  in  [44,  45,  47]  for  creep 
loads  of  smooth  specimens. 
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TABLE  V:  Cycle  number,  Characteristic  Size  and  Energy  Release  Rates,  A. 
due  to  Damage  Growth  for  two  Loading  Conditions. 


Configuration 

Cycle  Number 

Characteristic 

Size,  |im 

Energy  Release 
Rate,  J/m2 

Conditions  C1 : 

1 

27,000 

42 

2 

50,000 

52 

17.5 

3 

150,000 

80 

38.0 

4 

300,000 

120 

73.0 

Conditions  C2: 

1 

10,050 

58 

2 

19,000 

87 

32.0 

3 

30,000 

117 

105.0 
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Fig.  15  (a)  Illustration  of  a  displacement  response  at  point  x,  due  to  a 

discontinuity  at  %  given  by  expression  (5),  (b)  Schematic  showing  the 
approximation  of  discontinuities  within  a  rectangle  (see  text  for  details) 
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Fig.  16  Growth  rates  of  the  characteristic  size,  isotropic  expansion  and 

distortion  of  the  damage  zone  vs.  energy  reiease  rate  (Conditions  Cl ) 
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Fig.  17  Growth  rates  of  the  characteristic  size,  isotropic  expansion  and 

distortion  of  the  damage  zone  vs.  energy  release  rate  (Conditions  C2). 
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